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ABSTRACT 

Core collapse supernovae of massive (> 8 Mq) stars are formed primarily in OB associations and 
help blow giant superbubbles, where their collective shocks accelerate most of the Galactic cosmic rays. 
The spatial distribution of these stars is thus important to our understanding of the propagation of the 
observed cosmic rays. In order to better model the Galactic cosmic-ray distribution and propagation, 
we construct a three-dimensional spatial model of the massive-star distribution based primarily on 
the emission of the H II envelopes surrounding the giant superbubbles which are maintained by the 
ionizing radiation of the embedded O stars. The Galactic longitudinal distribution of the 205 /im N 
II radiation emitted by these H II envelopes is used to infer the spatial distribution of superbubbles. 
We find that the Galactic superbubble distribution is dominated by the contribution of massive star 
clusters residing in the spiral arms. 

Subject headings: cosmic rays: propagation — elementary particles — gamma rays: theory — inter- 
stellar medium - nucleosynthesis — supernovae 
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1. INTRODUCTION 

Supernovae (SNe) created by the core collapse of mas- 
sive stars have long been recognized (e.g., Baade & 
Zwicky 1934; Zwicky 1939; Hayakawa 1956; Ginzburg 
& Syrovatskii 1964; Zatsepin & Sokolskaya 2007) as the 
major source of Galactic cosmic-ray acceleration. How- 
ever, in the eighty years since then, subsequent obser- 
vations have demonstrated that collapse supernovae are 
quite complex. Presently, the principal categories of core- 
collapse supernovae (CCSNe), classified by their optical 
spectra, are: II- P (plateau), ILL (linear), Iln (narrow 
emission lines), lib (transitional), lb, and Ic (e.g., Fil- 
ippenko 1997). The minimum stellar mass necessary to 
produce a core collapse supernova has been found to be 
8 ± 1 Mq from the discoveries of red supergiant pro- 
genitors in Hubble Space Telescope (HST) pre-explosion 
images of SN II-P (e.g., Smartt 2009). At solar metal- 
licity an 8 Mq star corresponds to the spectral classifi- 
cation of B3 V (Drilling & Landolt 2000). We estimate 
the main-sequence age of such a 8 M star « 3.5xl0 7 
yr via interpolation from a tabulation of mass-dependent 
main-sequence ages (Massey & Meyer 2001). 

The most massive CCSNe progenitor most likely is as- 
sociated with SN Ic. The progenitors of SN Ic are very 
massive stars which have been stripped of both their 
hydrogen and helium envelopes (e.g., Filippenko 1997). 
However stellar progenitors have not been identified on 
SN Ic pre-explosion images (e.g., Smartt 2009). SN Ic 
models require the core collapse of either rotating sin- 
gle stars initially - 40-120 M (Georgy et al. 2009), 
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which have lost their stellar envelopes prior to collapse 
via stellar winds (e.g., Gaskell et al. 1986), or of ini- 
tially less massive (~ 25 M©) primaries (Dessert et al. 
2012) in close binary systems which have lost their stel- 
lar envelopes prior to collapse to smaller companion stars 
through Roche lobe overflow or common envelope evolu- 
tion (e.g., Pazyhski 1967; Podsiadlowski et al. 1992). At 
solar metallicity zero-aged 25 and 120 M Q stars corre- 
spond to spectral classifications of 08 V and 03 V; these 
are short lived with main-sequence ages of 6.4xl0 6 and 
2.6xl0 6 yr respectively (Massey 2003). 

These SN II-P and SN Ic constitute 48 ± 6% and 15 
± 4% respectively of CCSNe in spiral galaxies (Smith et 
al. 2011). 

These CCSNe progenitors are not distributed uni- 
formly in spiral galaxies. It has been known for close 
to a century that massive stars cluster (Kapteyn 1914). 
This is not surprising since the majority of stars are born 
in stellar clusters embedded in molecular clouds (McKee 
& Williams 1997; Lada & Lada 2003; Portegies Zwart et 
al. 2010). Moreover, the vast majority (70-80%) of all 
O stars reside still in the stellar clusters in which they 
were born (Gies 1987; Mason et al. 1998; Neugent & 
Massey 2011). The giant molecular clouds creating the 
massive CCSNe progenitor stars usually produce ~5 co- 
eval episodes of star formation separated by ~3.5x 10 6 yr 
(Blaauw 1964; McKee & Williams 1997). Together these 
sub-clusters constitute what is called an OB association 
(OB A) (Blaauw 1964). 

Analyses of Hipparcos proper motion measurements 
have found that the upper limit on internal OBA ve- 
locity dispersions is ~ 3 km s _1 (de Zeeuw et al. 1999). 
This modest velocity dispersion is the fundamental rea- 
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son why CCSNc occur primarily in supcrbubblcs (SB) 
since even an 8 M SN II-P stellar progenitor does not 
travel too far (~100 pc in 3.5xl0 7 yr) from its birth site 
before it dies as a CCSNe. Therefore the core-collapse 
deaths of (120-8 M©) progenitor stars occur within ~10 
to 100 pc of their birth sites. Each stellar death releases 
<~ 10 51 ergs of energy in the surrounding interstellar 
medium with a mean pre-supernova wind contribution of 
less than 10% (Leitherer at al. 1999). Consequently, such 
correlated supernovae create giant (> 150 pc) cavities of 
hot, tenuous plasma, or superbubbles, rather than many 
smaller isolated remnant bubbles (e.g., McCray & Snow 
1979; Tomisaka & Ikeuchi 1986; McCray & Kafatos 1987; 
Mac Low & McCray 1988; Tenorio-Tagle & Bodenhcimcr 
1988; Korpi et al. 1999; de Avillez & Breitschwerdt 2005; 
Stil et al. 2009). In such superbubbles cosmic rays are 
accelerated by CCSNe shock waves expanding through 
tenuous plasmas heavily enriched by supernova ejecta, 
synthesized by previous generations of CCSNe (Higdon 
et al. 1998; Higdon & Lingenfelter 2003, 2005; Ogliore et 
al. 2009; Rausch et al. 2009). 

A concrete example of a superbubble is the Orion- 
Eridanus superbubble. Heiles et al. (1999) term this su- 
perbubble the "Rosetta Stone of Superbubbles." Its over- 
all size is ~250 pc, it is filled with a x-ray (^2xl0 6 K) 
emitting plasma, and it is powered by stellar winds and 
CCSNe of massive stars created in earlier star-formation 
episodes in the Orion OBI association (e.g., Heiles et 
al. 1999; Heiles 2001). This large plasma volume is sur- 
rounded by a denser partial, expanding (^15 km s _1 ) 
shell which is ionized by the O stars of the Orion OBI 
association (Reynolds & Ogden 1979). This H II partial 
shell is in turn enveloped by an expanding H I shell of 
mass ~3xl0 5 M Q (Brown et al. 1995). 

In the present investigation we create a three- 
dimensional Galactic distribution of OB associations and 
their surrounding superbubbles. 

1.1. Outline 

Although the actual determination of the Galactic spa- 
tial distribution of superbubbles is straightforward, all 
the details involved make the derivation lengthly and 
complicated. 

In Section 2 we investigate the relationship between 
CCSNe occurrences and sites of massive star formation 
in other spiral galaxies. Of particular interest is the study 
by Anderson et al. (2012) of the correlation of SN II oc- 
currences and extragalactic star-formation sites as traced 
by near ultraviolet continuum emissions. Such emissions 
are sensitive indicators of the presence of moderate-mass 
(8 - 20 M ) SN II stellar progenitors. 

In Section 3 we discuss the spectral signatures spe- 
cific to superbubbles. We find that most appropriate SB 
signature turns out to be N II far infrared lines which 
have been resolved by the Far Infrared Absolute Spec- 
trophotometer (FIRAS) aboard the Cosmic Background 
Explorer COBE (Fixsen et al. 1999). We expect that 
such far infrared lines are radiated by ^10 2 pc sized (1 
- 10 cm~ 3 ) photoionized envelopes surrounding the hot, 
tenuous superbubbles. These H II envelopes absorb the 
majority of ionizing radiation emitted by the embedded 
OB associations (e.g., McKee & Williams 1997). 

In Section 4 we first consider a simple three- 



dimensional N II emissivity model and simulate the 
FIRAS 205 /im intensities, as a function of Galactic lon- 
gitude, assuming that the spatial distribution of Galactic 
superbubbles follows that of an axisymmetric disk popu- 
lation. This modeled FIRAS 205 /im intensity, displayed 
in Figure 1, poorly reproduces the measured FIRAS 
205 /xm intensity. More complicated source models are 
clearly required. 

Therefore in section 5 we present three-dimensional 
N II emissivity distributions based on models consist- 
ing of four spiral arms (e.g., Vallee 2008). We find that 
such spiral-arm distributions can reproduce the mea- 
sured (Fixsen et al. 1999) FIRAS 205 /an intensity once 
the contribution of two nearby (~ kpc) H II regions are 
included. 

In Section 6 we discuss the fact that the SB spatial 
distribution does not scale linearly with the N II 205 /im 
emissivity distribution due to the presence of a large- 
scale radial gradient in the ratio of interstellar nitrogen 
to hydrogen abundances. Here we present the spatial 
distribution of solely the massive star population. This 
distribution represents the best measure of the Galac- 
tic spatial distribution of superbubbles and their their 
embedded OB associations. 

In Section 7 we interpret the very bright 205 /jm 
FIRAS features at Galactic longitudes, around 80° and 
265°, not as a spiral-arm tangents, but as nearby H II 
sources. We model the strong 205 /im FIRAS feature 
at Galactic longitude, « 80° as a H II envelope of a su- 
perbubble, excited by a rich embedded OB association, 
Cyg OB2. We model the other 205 /im FIRAS feature at 
Galactic longitudes, w 265°, as a component of the Vela 
Nebula, a nearby (< 0.35 kpc) H II region, primarily 
excited by a single close O star, 7/-Puppis. 

Finally, in Section 8 we summarize the results of this 
investigation. 

2. THE RELATIONSHIP BETWEEN CCSNE 
OCCURRENCES AND STAR FORMATION SITES 

The majority (<~70%) of CCSNe are expected to be mid 
and early B stars, since massive stars follow a Salpetcr 
initial mass function (e.g., Bastian et al. 2010) and the 
minimum CCSNe progenitor mass is 8 M Q (e.g., Smartt 
2009). Yet Galactic OB A membership is invariably de- 
termined by bright O stars (Gies 1987; Mason et al. 1998; 
Mafz-Apellaniz et al. 2004). As noted by Maiz-Apellaniz 
et al. "the quality of the data regarding membership to 
a cluster or an OBA is extremely varied and for stars lo- 
cated at farther distances information is usually poorer." 
Cluster membership of the early and mid B stars, which 
constitute the majority of the CCSNe stellar progenitors 
are seldom investigated. In this section we discuss the 
evidence for stellar clustering of lower-mass, B-star (< 
20 M ) CCSNe progenitors. 

James & Anderson (2006) investigated the correlation 
of extragalactic CCSNe to star formation sites using Ha 
CCD images of the parent galaxies. In order to moder- 
ate contributions of faint star-formation sites they em- 
ployed all the pixels of an Ha image in their analysis. 
They sorted all the pixels of a digitized image in order 
of increasing pixel counts starting with the most nega- 
tive sky pixels (created by background subtractions) and 
ending with the pixel corresponding to the brightest star- 
formation region. The initial negative pixel values are 
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set to zero. Then a cumulative distribution is created 
from the ranked sequence and is normalized to the total 
Ha flux of the image. Thus a NCR (their acronym for 
the normalized cumulative rank pixel function) value of 
corresponds to background sky values, and NCR of 1 
corresponds to brightest pixel in an image. 

If the CCSNe progenitor was a member of the massive- 
star population which creates the tracer (here Ha) emis- 
sion, it is equally likely to come from any part of the 
NCR distribution (James & Anderson 2006). Thus over 
a long period of time the large number of pixels consti- 
tuting the 0.0 to 0.1 NCR segment contribute as much to 
total tracer emission as the few bright pixels comprising 
the 0.9 to 1.0 segment. In such a case the NCR distribu- 
tion will be uniform and have a mean value of 0.5 (James 
& Anderson 2006). 

Using the above NCR pixel statistic, Anderson et al. 
(2012) investigated the relationship between CCSNe and 
star-formation sites, traced by Ha emission, for a large 
sample of 260 (163 SN Type II and 96 Type Ib/c) CC- 
SNe. Their effective pixel size corresponds typically to a 
dimension of 130 pc. In view of their sample size they 
applied the NCR statistic to the six sub-type classifica- 
tions (II-P, II-L, Iln, IIP, lb, & Ic) applicable to CCSNe. 
Their results are listed in Table 1, where n refers to the 
number of SNc in the sample. Using the Kolmogorov- 
Smirnov (KS) test they found that three NCR distribu- 
tions, SN lib and SN II-L, as well as SN Ic, agree well 
with a massive-star population as traced by Ha emis- 
sion. Yet, the KS test demonstrates that two Type II 
sub-types, SN II-P (the most common CCSNe) and SN 
Iln, as well as SN lb, are not consistent with a massive- 
star population as traced by Ha line emission. However, 
as we shall see, Ha emission is not the sole tracer of 
massive stars. 



Table 1. Anderson et al. (2012) 



Correlations of CCSNe with Ha Tracer 


SNc Class 


n 


<NCR> Ha 


lb 


39 


0.32 ± 0.040 


Ic 


52 


0.46 ± 0.040 


II-L 


13 


0.37 ± 0.100 


lib 


13 


0.40 ± 0.100 


II-P 


58 


0.26 ± 0.039 


Iln 


19 


0.21 ± 0.065 



Anderson et al. also investigated further the correla- 
tion of CCSNe and star formation sites, but now used 
near ultraviolet (NUV 1750-2750 A) images taken by 
the Galaxy Evolution Explorer (GALEX) (Martin et 
al. 2005). NUV continuum emission is particularly well 
suited to identify > 10 Myr old stellar clusters, which ra- 
diate negligible Lyman continuum (Lyc) emission. Such 
UV emission provides a direct measure of star formation 
on timescales ~10 8 yr (Martin et al. 2005). Note that a 
~10 8 yr lifetime corresponds to 5 M B5 V star (Massey 
& Meyer 2001). Using NUV emissions they found for 
the SN II-P (50 SNe) <NCR> NUV = 0.50 ± 0.045 and 
SN Iln (18 SNe) <NCR>aw = 0.38 ± 0.065. The KS 
test shows that these two NCR distributions are consis- 
tent with that expected for flat distributions, and, con- 
sequently, these two SN II sub-types track well stellar 



clusters containing older (10-100 Myr), but still massive 
(5-20 Mq) stars. 

Thus analyses of sites of extragalactic CCSNe demon- 
strate that the great majority of these CCSNe occur in 
or close to (< 150 pc) of stellar clusters of massive stars 
all with the potential to generate superbubbles. 

3. MODELING OF THE GALACTIC SUPERBUBBLE 
SPATIAL DISTRIBUTION 

We are interested in radiation signatures specific to 
superbubbles. The clearest indicator of superbubbles is 
soft x-ray emission from large (> 10 2 pc) cavities filled 
with hot (^10 6 K) plasmas. Yet photoelectric absorp- 
tion by intervening interstellar gas hinders the detection 
of such soft x-ray emissions beyond the solar vicinity. 
In their seminal study of Galactic OB associations Mc- 
Kee & Williams (1997) noted that moderate-density (~ 3 
cm~ 3 ) photoionized H II envelopes surrounding hot, ten- 
uous superbubbles cores absorb the majority of the ion- 
izing radiation emitted by embedded association O stars. 
Therefore the determination of the spatial distribution of 
these H II envelopes provides an accurate measure of the 
distribution of OBA, the primary Galactic CCSNe sites. 

Previously Anantharamaiah (1986) proposed such H 
II envelopes as the primary Galactic source of the dif- 
fuse low frequency (~325 MHz) radio recombination lines 
(RRL). Subsequent analyses of a larger RRL data set by 
Roshi & Anantharamaiah (2001) strengthened the inter- 
pretation of moderate-density H II envelopes as the origin 
of low-frequency RRL: first, RRL emissions exhibit large 
fluctuations as a functions of longitude, which is expected 
as lines of sight intercept discrete H II envelopes of dif- 
fering properties, second, the path lengths through the 
individual clumps range from 20 to 200 pc, and, third, 
the mean electron densities, n e , are ~1-10 cm~ 3 . Con- 
sequently, for the reasons presented here we will employ 
observations of H II envelopes as tracers of Galactic su- 
perbubbles. 

In their investigation of OB associations McKee and 
Williams (1997) employed far infrared N II line emis- 
sions together with thermal radio continuum radiation 
to determine the Galactic Lyc emission of OB associa- 
tions. We will use far infrared N II line emissions to 
determine the spatial distribution of H II envelopes of 
superbubbles. Recently, Steiman-Cameron et al. (2010) 
have employed such N II line emissions to determine the 
Galactic spatial distribution of photoionized gas. Our 
study differs from this investigation in several important 
respects, which will be discussed later. 

The rationale for using far infrared N II line emissions 
to trace H II emission is the following. Since the ioniza- 
tion potentials of H I and N I are 13.6 eV and 14.5 cV, 
H II and N II regions are essentially coincident. Thus we 
use the N II far infrared lines at 122 /zm & 205 yum as 
observed by the Far Infrared Absolute Spectrophotome- 
ter (FIRAS) aboard the Cosmic Background Explorer 
COBE (Fixsen et al. 1999) to infer the spatial distri- 
bution of these moderate-density H II envelopes. Note 
that the ratio of 122 /im line flux to 205 /im line flux is 
density dependent and it rises with increasing electron 
density, n e . In the limit of very low densities, where all 
the nitrogen ions are in the ground state, the ratio of 
the 122 ^m line flux to 205 /im line flux is 0.7; the ratio 
rises from 1 to 1.6 for n e between 10 and 35 cm~ 3 for 
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uniform-density H II regions (Wright ct al. 1991). Yet 
at n e of 10 2 cm~ 3 the ratio increases to 3 and at even 
higher densities it saturates at 10 (Rubin 1985). How- 
ever Fixsen et al. found that FIRAS instrumental noise 
precluded a detailed derivation of 122 pm line flux but 
the overall Galactic disk mean ratio of 122 pm line flux 
to 205 pm line flux is 1.1 ± 0.1. Thus the FIRAS N II 
flux ratio infers that typical electron densities similar to 
the values expected (e.g., McKee & Williams 1997; Roshi 
& Anantharamaiah 2001) for H II envelopes. 

Where are such superbubbles expected to be located 
in the Galaxy? In their investigation McKee & Williams 
(1997) made the simplifying assumption that the surface 
density of OB associations could be approximated by an 
axisymmetric Galactic disk population, S(p) oc e~ p / Hp , 
where S represents the population's surface density and p 
a source planar distance from the Galactic center. They 
determined that H p = 3.34j^'24 kpc when they used an 
investigation of luminous Galactic H II regions (Smith, 
Biermann, & Mezger 1978), which found that all bright 
disk H II regions are restricted to p m in < P < Pmax where 
Pmin — 0.39R S , Pmax = 1-3R S , and R s is the distance of 
the Sun from the Galactic center. We now model the 
longitudinal FIRAS 205 p intensity distribution corre- 
sponding to an axisymmetric truncated exponential disk. 

4. CALCULATION OF THE SPATIAL DISTRIBUTION OF 
THE N II EMISSIVITY FOR A TRUNCATED 
AXISYMMETRIC DISK SOURCE MODEL 

Such a source distribution is best represented in a 
Galactocentric Cartesian coordinate system, where the 
y-axis is aligned along the Sun-Galactic Center line, and 
the x-axis, which passes through the Galactic Center 
(GC), is parallel to the galactic longitude 1 = 90°. There 
a Galactocentric cylindrical {p, 9, z) coordinate system 
is parameterized as p = ^/(x 2 +y 2 ), 9 = tan _1 (y/x), and 
z is measured normal to the plane. We employ the above 
S(p) to represent the disk N II 205 ^m emissivity and 
then determine the corresponding longitudinal distribu- 
tion of the N II 205 pm intensity. The N II intensity 
must be calculated in Earth-centered coordinates (r, 1, 
b), where r is distance, 1 is Galactic longitude, and b 
is Galactic latitude, which are related to the previously 
defined Galactocentric cylindrical coordinates (p, 9, z), 



p = \J r 2 cos 2 b + R 2 — 2R s cosbcosl, 

_i ,R S — rcosbcosl 

tan ( — ), z = rsinb. 

rcosbsinl 



(1) 



The N II 205 pm line emissivity, e D , in units of erg s 1 
cm~ 3 , is calculated from S(p) in the following way, 



\p) = L NII 2Q5p,mS(p)P z (z)/A 



D, 



(2) 



where the Galactocentric planar distance p(r, 1, b), is a 
function of Earth-centered (r, 1, b), Ljyj 7205^1 is the total 
Galactic 205 pm luminosity in erg s , P z (z) represents 
the height dependence of the H II SB envelopes, and 
represents the source-weighted Galactic-disk area. Once 
P z (z) is defined, the disk emissivity can be used to cal- 
culate a line energy flux, F, 



Jlmin Jb m i n JO 



,D 



(p(T,l,b)) 
Aitr 2 



r 2 drdlcosbdb, (3) 



in units of erg cm -2 s . 

Using COBE, Bennett et al. (1994) performed a far- 
infrared spectral line survey from over 99% of the sky 
with a 7° beam, resolving the Galactic longitudinal dis- 
tributions of 122 and 205 pm N II lines, see their Fig. 
2d and 2e. Subsequently, Fixsen et al. (1999) analyzed 
further the FIRAS data and produced a factor of 2 im- 
provement in the signal to noise ratio of the line emis- 
sions as well as a more complete sky map of the N II 205 
pm line emissions. Consequently, we employ the later 
analyses of Fixsen et al. (1999) in comparisons with our 
intensity models. In their Fig. 5e Fixsen et al. (1999) 
plot the mean intensity in units of W m~ 2 sr _1 of the N 
II fine-structure ground-state transition at 205 pm aver- 
aged over Galactic longitudes (centered on the Galactic 
plane) of Al = 5° and Ab = 1°. As noted by Fixsen et al. 
(1999) the 7° FIRAS beam does not resolve the latitu- 
dinal distributions of the far infrared lines, consequently 
we compare the FIRAS N II 205 pm line measurements 
with our models of the mean line intensity, averaged over 
Ab =1° and Al = 5° in the following way, 



d 2 F 
dldb 



1 

: Ab 



-3.5 



-3.5° 



\p(r,l,b)) 
4irr 2 



r 2 drcosbdb, (4) 



in units of erg cm' 

Before we calculate our modeled disk intensity, we 
must address how we represent P z (z). Here we assume a 
Gaussian distribution as a function of z, 



2-K(J z 



(5) 



As noted above, FIRAS observations did not resolve 
the latitudinal extents of the far infrared line fluxes, 
yet the low-frequency RRL observations (Roshi & Anan- 
tharamaiah 2001) found that in the inner Galaxy the 
low-frequency RRL intensities as a function of Galactic 
latitude, b, possess full-width half-maxima (FWHM) of 
~1.8°. In order to reproduce such FWHM transverse the 
Galactic plane our modeling of the diffuse Galactic N II 
205 pm line emissions required a z w 0.15 kpc, signifi- 
cantly greater than <~ the 0.045 kpc, scale height of OB 
stars (e.g., Reed 2000) or the 0.035 kpc, scale height of 
molecular clouds (e.g., Stark & Lee 2005). But our <r z 
of 0.15 kpc is consistent with the scale height expected 
(McKee & Williams 1997) for H II envelopes. 

Our Figure 1 shows a modeled FIRAS 205 pm mean 
intensity distribution corresponding to an axisymmetric 
truncated exponential disk with a radial scale, H p , of 2.5 
kpc. Note that the radial gradient of the N II emissivity 
distribution is steepened (e.g., McKee & Williams 1997) 
by the presence of a radial gradient in the interstellar 
nitrogen to hydrogen abundance ratio. Thus a massive 
star distribution with a radial scale, R p of 3.4 kpc, creates 
a steeper N II emissivity distribution with an H p of 2.5 
kpc due to the existence of a large-scale Galactic gradient 

in N/H, c - p / H p /H with a Rp /H of 10 kpc (Daflon & 
Cunha 2004). We normalized this modeled FIRAS 205 
pm intensity to the measured line intensity at 1 = 30° 
which required a total Galactic N II 205 pra luminosity of 
1.85xl0 40 erg s" 1 . Inspection of this Figure shows that 
such a modeled intensity, created by an axisymmetric 
disk distribution, decreasing monotonically in p, falls off 
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too slowly with Galactic longitude > 1 30° | and fails to 
reproduce the observed FIRAS longitudinal intensities. 
More complicated emissivity models are necessary. 

5. CALCULATION OF THE SPATIAL DISTRIBUTION OF 
THE N II EMISSIVITY FOR SPIRAL- ARM SOURCE 
MODELS 

Based on optical observations Morgan et al. (1952) 
were the first to discern that H II regions traced out spiral 
arms in the Milky Way Galaxy. In a later groundbreak- 
ing study, using both Ha and high-frequency RRL ob- 
servations of H II regions, Georgelin & Georgelin (1976) 
found that the Galactic H II distribution delineated a 
spiral structure comparable to the large-scale spiral con- 
figurations observed in other galaxies and the majority 
(80%) of the luminous Galactic H II regions fell along 
two symmetrical pairs of arms with pitch angles of 12°. 
Subsequent studies of spiral galaxies have demonstrated 
that the principal spiral-arm tracers are massive star for- 
mation sites defined by massive stars, photoionized gas 
and giant molecular clouds (e.g., Baade 1963; Boulanger 
et al. 1981; Rumstay & Kaufman 1983; Russeil 2003; Hou 
et al. 2009; Efremov 2011). 

Since that first pioneering investigation by Morgan et 
al. (1952) Galactic spiral arms have been traced primar- 
ily by H II regions (e.g., Georgelin & Georgelin 1976; 
Downes et al. 1980; Caswell & Haynes 1987; Watson et 
al. 2003; Russeil 2004, 2007; Paladini et al. 2004; Hou et 
al. 2009) as well as by giant molecular clouds (e.g., Co- 
hen et al. 1986; Dame et al. 1986; Grabelsky et al. 1988; 
Digel et al. 1990; May et al. 1997; Heyer et al. 2001; 
Dobbs et al. 2006; Stark & Lee 2006; Nakanishi & So- 
fue 2006; Hou et al. 2009). However, other tracers such 
as H I (e.g., Englmaier & Gerhard 1999; Nakanishi & 
Sofue 2003; Levine et al. 2006) as well as infrared emis- 
sions (e.g., Hayakawa et al. 1981; Bloemen et al. 1990; 
Drimmel & Spergel 2001; Benjamin et al. 2005) also have 
been used to determine some spiral-arm properties, par- 
ticularly, the locations of spiral-arm tangents. 

In a detailed study evaluating and critically weigh- 
ing the past evidence of Galactic spiral structure, Vallee 
(2008) has compiled key spiral-arm parameters and their 
uncertainties and, consequently, he determined up-to- 
date locations of four logarithmic spiral arms. We em- 
ployed Vallee's spiral-arm parameters in our models of 
N II line FIRAS intensities. We have expanded Vallee's 
model by the inclusion of the fall offs of the spatial den- 
sity of H II envelopes parallel to the spiral arms and nor- 
mal to them. In the present section we now model the 
spatial distribution of the superbubble spiral-arm popu- 
lation. 

Vallee's spiral-arm model can be described in the 
Galactocentric cylindrical coordinate system, presented 
earlier, by 

p{9) = p min e k{e - e ^ (6) 

Pmin — P — Praax ) 

where k = tanp, and p is the pitch angle. From 
Vallee's Fig. 1 we determined the following arm param- 
eters for the starting points of the arms: p m in of 2.9 
kpc, o of 70° (Norma-Cygnus arm), 160° (Perseus arm), 
250° (Sagittarius-Carina arm), and 340° and (Scutum- 
Crux arm). For the spiral-arm pitch angles, p, we em- 
ployed a value of 13.5°, the unweighted mean value from 



Vallee's Table 1 for three spiral arms (Norma-Cygnus, 
Sagittarius-Carina, and Perseus). In order to reproduc43 
the observed FIRAS 205 fim line flux at 1 = 310° we 
used a somewhat greater pitch angle, p, of 15.5° for the 
Scutum-Crux arm in order to place its tangent at this 
longitude. We employ Vallee's p m ax of 35 kpc. Our ren- 
dition of Vallees arm model is shown in our Figure 2. 
We use R s = 7.6 kpc, the distance of the Sun from the 
Galactic center. 

In our investigation we employ the above Eq. (6) to 
represent the medians of the spiral arms and introduce 
three parameters to construct a three-dimensional spa- 
tial distribution for a spiral-arm population. First, we 
modeled the decrease of the arm population with Galac- 
tocentric radius with an exponential distribution 

P p {p) = e-"I H ". (7) 

McKee & Williams (1997) noted that such an exponen- 
tial distribution with H p « 3.3 kpc reproduced well the 
axisymmetric Galactic spatial distribution of young stel- 
lar populations as well as that of giant molecular clouds. 
Second, following Taylor and Cordes (1993) modeling of 
the spiral-arm component of Galactic thermal electrons, 
we represent the fall off of arm populations across an arm 
with distance, A, transverse to the arm medians of the 
modeled arms, as Gaussian distributions, 

P A (A) = e -°- 5A2 /^ ( 8 ) 

with a fixed dispersion, a a ■ Third, we employed a Gaus- 
sian distribution, Eq. 5, as well as a z of 0.15 kpc, to 
represent the fall off of H II envelopes with height, z, 
transverse the Galactic plane. Now we derive realistic 
values of H p and a a from model comparisons with the 
Galactic N II 205 pm longitudinal intensities of Fixsen 
et al. (1999). 

For every pair of chosen arm parameters, H p and a a, 
the corresponding planar arm distribution is calculated 
using the following approach. In the Galactocentric co- 
ordinate system, defined earlier, coordinate pairs (pk, 0k) 
are calculated along an arm median from Eq. 1. Then 
at each pk, Eq. 7 is used to calculate relative arm source 
density, P p (pk). At these latter locations, further posi- 
tions, (pj,8j), are calculated at fixed distances, Aj, pro- 
jected normal to each arm. At these positions, Eq. 8 is 
used to calculate relative arm density, Pa(Aj), at each 
Aj ranging non-uniformly from 0.01 kpc to 5 kpc. At the 
start (pfe = Pmin) and the termination (p^ = p ma x) of an 
arm, where our knowledge of arm properties is particu- 
larly fragmentary, each Pa(A^) is projected in a circular 
arc, and a function of Aj, as is illustrated in our Fig. 
2 forAj = 0.5 kpc. Thus the relative arm source densi- 
ties, P p (Afe)xPA(Aj), are calculated for a sample of 10 5 
locations in the vicinities of the arms. 

The fractional contribution of each spiral arm to the to- 
tal N II Galactic luminosity ( in erg s _1 ) is found by inte- 
grating the product, P p (Afc)xPA(A :/ ), over the Galactic- 
disk area and normalizing the surface integrals to N II 
luminosity fractions discussed below fx, where X repre- 
sents NC (Norma-Cygnus), SA (Sagittarius-Carina), SC 

1 Using a pitch angle, p. of 13.5° for the Scutum-Crux spiral 
arm produced arm tangents at 1 = 26° and 1 = 315°; changing p 
to 15.6 circ placed arm tangents at 1 = 27° and 1 = 311°. 
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(Scutum-Crux), and P (Perseus). In order to use such 
modeled arm distributions, interpolants, Sx(p, #), were 
created (for each arm X) to permit the calculation of 
planar source densities as functions of arbitrary Galac- 
tocentric coordinates (p, 6), 

Sx(p,0) = f x P p (p k (p,6)) x P A (p 3 ,0 3 )/A x , 

A x = / P p (p k (p,6))P±(Aj{p,0))pdpd0. 
Jo Jo 

Note that each Sx{p, 0) distribution, when integrated 
over the disk area normalizes to ix, and, thus, the sum 
of these four integrated distributions equals unity. 

We now employ these Sx {p, 0) to represent the spiral- 
arm contribution of the N II 205 pm emissivity and then 
determine the corresponding longitudinal distribution of 
the N II 205 pm intensity. The N II 205 pm line emissiv- 
ity in each arm, e x in units of erg s _1 cm~ 3 , is calculated 
from Sx {p, 0) in the following way, 

e x ( Pl 9, z) = L NII205fim S x (p(r, I, b),0(r, I, b))P x (z), (9) 

where the Galactocentric planar distance, p(r, 1, b), and 
the polar angle, 0(r, 1, b), are functions of Earth-centered 
(r, 1, b), Ljv/7205Ai"i i s the total Galactic 205 pm lumi- 
nosity in erg s" 1 . Using Eq. 5 to model P z (z), these arm 
emissivities can be used to calculate mean line intensity, 
averaged over Ab = 1°, in the following way, 

w=W_, 5 c I ~^^ r drcosbdb > (10) 

in units of erg cm~ 2 s _1 sr _1 . 

In Fig. 3 intensities, d 2 F/dldb, summed over all four 
arms, calculated from Eq. 10, are shown for a total Galac- 
tic N II 205 pm luminosity of 1.4x 10 40 erg s _1 and nom- 
inal spiral arm parameters of a a = 0.5 kpc and a H p 
= 2.4 kpc for four sets of relative arm luminosities, Isa, 
fsc, f/vc, and fp. The calculations illustrated in Fig. 3 
demonstrate that the equal spiral-arm luminosities can- 
not reproduce the observed FIR AS 205 pm longitudinal 
intensity fall off at 1 > 35° and 1 > 335° in the N II 205 
pm intensities and, consequently, the contributions of the 
Norma-Cygnus and Sagittarius-Carina arms must be de- 
creased by a factor of <~ 2 relative to the contributions 
of the Scutum-Crux, and Perseus arms. In their spiral- 
arm model of far infrared N II line emissions Steiman- 
Cameron et al. (2010) noted similar behavior. 

Such differences in the arm contributions have also 
been identified in analyses of other spiral-arm trac- 
ers. Drimmel and Spergel (2001) modeled near in- 
frared and far infrared continuum emissions measured 
by DIRBE/COBE and they found that two spiral 
arms (Scutum-Crux and Perseus) dominate the infrared 
spiral-arm signatures and that the contribution of the 
Sagittarius-Carina arm must be reduced relative to the 
contribution of the Scutum-Crux arm. Such a two- 
armed continuum infrared spiral signature is supported 
by a recent derivation (Benjamin et al. 2005; Ben- 
jamin 2008) of the stellar source distribution from the 
Spitzer / GLIMPSE mid-infrared survey of the Galactic 
plane. Recently Efremov (2011) noted that both the 
Sagittarius-Carina and Norma-Cygnus arms are better 



defined in H I than the other two arms but young stellar 
populations exist "beyond a doubt" in the Sagittarius- 
Carina arm. Moreover such arm differences have been 
identified in other spiral galaxies. For example, the 
derivation of spatial distributions of H II regions in the 
two-armed Sc galaxy, NGC 628, by Kennicutt & Hodge 
(1976) showed that the mean H II densities differ by 1.8 
between the two arms. 

In view of this discussion and the our modeled inten- 
sities illustrated in Fig. 3, in the remainder of this study 
we employ: is A — 0.18, Lvc = 0.18, fsc = 0.28, and fp 
= 0.36. 

In Fig. 4 intensities, d 2 F/dldb, summed over all four 
arms, calculated from Eq. 10, are shown for a total Galac- 
tic N II 205 pm luminosity of 1.4x 10 40 erg s _1 for a fixed 
u A = 0.5 kpc and a range of radial arm parameters H p 
from 1.8 to 3.0 kpc showing that the value of 2.4 kpc gives 
the best fit. In Fig. 5 intensities, d 2 F/dldb, summed over 
all four arms, calculated from Eq. 10, are illustrated for a 
fixed total Galactic N II 205 pm luminosity of 1.4xl0 40 
erg s^ 1 for a fixed radial arm parameter, H p , of 2.4 kpc 
for a range of arm-width parameters, a a-, from 0.25 to 
0.75 kpc showing that the value of 0.5 kpc gives the best 
fit. Also we include in these last two figures our estimates 
of the contributions of the two nearby H II regions: a 
~6°-sized Cygnus X diffuse H II region with its center at 
1 = 80° & b = 0°) as well as from 10°-sized Vela-complex 
H II region with its center at 1 = 262° & b = -2°). The 
properties of these two nearby H II regions are addressed 
in section 7. 

From these two figures, we find that H p = 2.4 ± 0.3 
kpc, and a a = 0.5 ± 0.1 kpc reproduce well the diffuse 
FIR AS ljNH205fmi longitudinal intensity, when the two 
strong local H II sources are included. These parameters 
this define our three-dimensional spiral-arm luminosity 
model for N II 205 p m emission, as shown in Fig. 6 to- 
gether with the individual contributions from each of the 
spiral arms, for a total spiral arm N II 205 pm luminosity 
of w (1.4 ± 0.15)xl0 40 erg s" 1 . 

Using these values and our nominal Galactic luminos- 
ity N II 205 pm of 1.4 xlO 40 erg s _1 our model results 
suggest that total Galactic 205 pm line luminosity in- 
side the solar circle (p < R s ) is 1.2xl0 40 erg s -1 . From 
their analyses of FIR AS line measurements Bennett et 
al. (1994) estimate, with a ~ 50% uncertainty, a some- 
what greater total 205 pm Galactic line luminosity inside 
the solar circle of 2.6xl0 40 erg s _1 . 

As can be seen in Fig. 4 and 5 all our models system- 
atically overproduce the observed N II 205 pm intensity 
along the tangent to the Sagittarius spiral arm at 55° < 
1 < 30°, where it is known (e.g., Georgelin & Georgelin 
1976; Forbes 1983) that a 5-kpc portion of the Sagittarius 
is "devoid" of luminous H II regions. 

From these figures we also see that the introduction of 
a galactic interarm population is not required. A lack of 
evidence for interarm superbubble population is perhaps 
surprising. Recent analyses of high-resolution Ha obser- 
vations of two nearby spiral galaxies shed some light on 
this thorny issue. Employing HST images, Lee, Hwang, 
& Lee (2011) cataloged 19,598 H II regions at a resolu- 
tion of 0.1" (5 pc) in M51, a system of two interacting 
galaxies, NGC 5194, a two-armed, grand-design spiral, 
and a small barred lenticular galaxy NGC 5195. In NGC 
5194 they determined the numbers of H II regions in 
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the arms, nucleus, and the interarm region, respectively 
12245, 4422, and 2636. They found that the relatively 
less luminous interarm H II regions contribute only 8% 
of total Galactic Ha H II luminosity of NGC 5194. Re- 
cently, Azimlu, Marciniak, & Barmby (2011) cataloged 
3961 H II in the neighboring Sb galaxy, M31, with a res- 
olution of 1" (3.8 pc). They determined that, although 
^40% of the H II regions occurred in the interarm region, 
the much less luminous interarm H II regions contributed 
only 15% of the total Galactic Ha luminosity. Based on 
such Ha observations we suggest that a modest (~10%) 
interarm SB contribution is difficult to discern due to 
uncertainties in the spiral-arm models as well as in the 
contribution of nearby H II regions. 

Recently Steiman- Cameron et al. (2010) published a 
spiral-arm source model to characterize COBE / FIRAS 
far infrared emissions for the C II 158 /im line as well 
as for the N II 205 (im line. Our approach differs from 
their study in several important respects. We focus our 
investigation solely on 205 /tm radiation from H II inter- 
stellar phases, assuming an H II envelope superbubblc 
origin, and, consequently, we used as ancillary data low- 
frequency RRL observations from Roshi & Ananthara- 
maiah (2001). This led to our choice of a z of 150 pc 
for the Gaussian dispersion for H II envelopes transverse 
the Galactic plane. Steiman-Cameron et al. did not em- 
ploy the H II envelope superbubble origin for N II gas. 
Instead they used a er z of 70 pc that is appropriate to 
the C II 158 /im latitude distribution which has been re- 
solved in a balloon experiment (Nakagawa et al. 1998). 
A significant component of the C II 158 /im line emis- 
sion may be emitted (e.g., Bennett et al. 1994; Roshi 
& Anantharamaiah 2001) by dense photodissociation re- 
gions of molecular clouds which possess a relatively small 
scale height transverse the Galactic plane (e.g. Stark & 
Lee 2003). Thus we chose a greater value for a z . Fur- 
ther, they convolved their modeled intensities over the 
FIRAS beam while we employed the published (Fixsen et 
al. 1999) area-averaged intensities. We also represented 
the geometry of the spiral-arm terminators differently. 

Nonetheless the COBE /FIRAS 205 fim data is robust. 
Although there exist significant differences in approach, 
the fundamental results of both investigations are the 
same: first, four- armed spiral models with differing rela- 
tive arm luminosities fit the FIRAS 205 /im longitudinal 
distributions, and, second, an interarm source contribu- 
tion is not required. 

6. CALCULATION OF THE SPATIAL DISTRIBUTION OF 
LYMAN CONTINUUM EMISSION 

Based on this model for the N II 205 itm emission, we 
now calculate the spiral arm density distribution of the 
Lyc emission from the O stars which produced it. As 
discussed above the Galactic Lyc does not scale simply 
with N II 205 /tm emissivity distribution due to the pres- 
ence (e.g., Shaver et al 1983) of a gradient, as a function 
of Galactocentric radius, in the interstellar abundance 
of nitrogen. Subsequently, Daflon & Cunha (2004) ana- 
lyzed the emissions of OB stars located at 4.7 kpc < p < 
13.2 kpc and derived a nitrogen abundance gradient of - 

0.042 dex per kpc or an equivalent of 10 kpc. Thus 

from the spiral-arm N II emissivity distributions H p = 
2.4 ± 0.3 kpc that reproduced the FIRAS 205 /im line in- 
tensities, we expect that 1/H^ /J w l/Rp /H + l/TA L p YC , 



and, consequently, the actual fall off in the Galactic Lyc 
emissivity, Hp YC = 3.2 ± 0.5 kpc. This is in excellent 
agreement with the axisymmetric source model of McKce 
and Williams (1997) where H p = 3.3 kpc. From a radial 

scale of Yl^ H of 10 kpc for the interstellar N to H ratio 
and nominal spiral- arm parameters (H p = 2.4 kpc and 
a A = 0.5 kpc) we find that the Galactic source weighted 
N/H = 1.78(N/H) . 

The resulting spiral-arm Lyc luminosity density model 
for Hp = 2.4 kpc and <ja = 0.5 kpc, normalized to a 
total Galactic Lyc luminosity of unity is shown in Fig. 7. 
This figure represents the best measure of the Galactic 
spatial distribution of superbubbles and their embedded 
OB associations. 

7. NEARBY SOURCES OF DIFFUSE 205 fjM N II 
EMISSIONS 

Early investigations (e.g., Weaver 1970; Jackson et al. 
1979; Vazquez et al. 2008) of young stars in the solar 
vicinity (< kpc) reported that the Sun is "located in an 
arm-like structure" terminated by the Cygnus region at 
one end at 1 ~ 80° and Vela at the other end at 1 ~ 264°; 
even then such a feature was termed (e.g., Weaver 1970) a 
"spur" or "offshoot" of the Sagittarius arm. Vallee refers 
to the structure as the "Orion armlet" or "Orion bridge" 
describing it as a "local aggregate of stellar clusters". 
This local structure is not viewed as a separate compo- 
nent of the global spiral-arm structure (e.g., Georgelin & 
Georgelin 1976; Taylor & Cordes 1993; Vallee 2008). Yet 
inspection of the observed FIRAS 205 /mi longitudinal 
intensity (Fixsen et al. 1999), illustrated in our Figs. 1, 
3, 4, 5, and 6, show that significant emission at 1 <~ 80° 
and 1 ~ 265° from this "Orion Spur" , or "local aggregate 
of stellar clusters" . 

Therefore in addition to the contribution of the large- 
scale, Galactic spiral arm distribution in this investiga- 
tion we consider the emission from these two nearby (~ 
kpc) H II regions, the Cygnus X region at 1 w 80° and the 
Vela Complex at 1 w 265°. A comparison of measured 
(Fixsen et al. 1999) FIRAS intensities at N II 205 /im 
(Fig. 5e) and at 122 /tm (Fig. 5f) show that the ratios 
of these two line intensities in the vicinities of 1 w 80° 
and 1 w 265° agree with the overall Galactic ratio. Con- 
sequently, we assume that the photoionized gas which 
radiates these lines has approximately the same electron 
density, ^3 cm' 3 , as the H II envelopes in the inner 
Galaxy which radiate the bulk of the spiral-arm 205 /tm 
N II line emission. 

7.1. Cygnus Region 

We interpret the bright 205 /tm FIRAS feature at 1 
w 80°, not as a spiral-arm tangent, but as a H II en- 
velope excited by a rich, nearby OB association, Cyg 
OB2. The Cygnus-OB2 stellar cluster is a visually ob- 
scured, but nearby, OB association at a distance of 1.45 
kpc (Hanson 2003) with its center at 1 = 80.1° and b 
= +0.9° (Humphreys 1978). In a groundbreaking study 
employing infrared observations in the J, H, and K bands 
Knodlseder (2000) estimated that this rich stellar cluster 
contained as many as 100 O stars and, clearly, it is the 
closest massive Galactic star formation site. Recently, 
employing Chandra x-ray point-source data as well as 
complementary optical and near IR photometry, Wright 
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et al. (2010) found that Cygnus OB2 encompasses an es- 
timated 75 O stars as well as ^1200 OB stars. There is a 
large spread of stellar ages, 2-10 Myr, in this rich stellar 
cluster (e.g., Hanson 2003; Drew et al. 2008; Wright al. 
2010; Comeron & Pasquali 2012). Most recently, from 
this massive star-formation site Fermi has detected and 
resolved a large 4° (~ 10 2 pc) superbubble, illuminated 
in 0.1 to 10 2 GeV 7-rays emitted by freshly-accelerated 
cosmic rays in the Cygnus Superbubble (Ackermann et 
al. 2012). This exciting new observation provides direct 
evidence that cosmic rays are accelerated in superbub- 
bles. 

Lozinskaya et al. (2002) investigated the kinematics of 
photoionized gas in Cygnus region via Ha measurements. 
In this study they posed a fundamental question - "where 
are the observational manifestations of the action of in- 
tense ionizing radiation from Cyg OB2 on the ambient 
gas?" Lozinskaya et al. propose that Cygnus X is the 
"observational manifestation" of Cyg OB2. Cygnus X 
is an extended H II region delineated by an irregular 
ring of optical Ha filaments ~ 8° in radius (Ikhsanov 
1961). Other investigates have also suggested Cygnus X 
and Cygnus OB2 are related (e.g., Ikhsanov, 1961; Veron 
1965). 

We propose that the H II region, Cygnus X, ionized by 
O stars in the embedded Cyg OB2, radiates the bulk of 
the 205 /im line emission resolved by FIRAS from 1 rj 
80°. We estimate the corresponding ionizing luminosity, 
Qo, for this association to be s» l.lxlO 51 s _1 , from 75 
young O stars, employing a time-averaged, ionized stellar 
luminosity relation^ and a initial stellar mass function 
with a slope of -1.35 (Salpeter 1955). The Cyg OB2 
association is very luminous, exceeding by a factor of 
~40, the ionizing luminosity, Qo, of the familiar Orion 
OB association of only 2.7x 10 49 s _1 (Williams & McKee 
1997). 

Lozinskaya et al. (2002) pointed out that the absence of 
intense Ha emissions, expected for such a large number 
of O stars, could be explained by strong visual absorp- 
tion, but corresponding thermal radio emissions must be 
observed regardless of visual obscuration. However, they 
remark that in the Cygnus X region radio sources are 
detected against a strong diffuse radio source ~ 6° in 
size centered on the Cyg OB2 association (e.g., Wend- 
ker 1970; Huchtmeier & Wendker 1977; Lozinskaya et 
al. 2002), In a 2.7 GHz radio survey an angular resolu- 
tion of 11 Wendker (1970) measured a total flux density, 
S v , of 5390 Jy for the total emission including resolved 
sources and solely 2260 Jy for the diffuse component of 
the Cygnus X. 

The ionizing luminosity, Qo, required to produce such 
diffuse thermal radio emission can be estimated (e.g., 
Afflerbach, Churchwell, & Werner 1997), 

7.6 xlO 48 S v v r 
Q ° ~ Ti' 2 { Ty ){ GWz } Wc> 

in units of s _1 where v is the radio frequency in GHz, T e 
represents the electron temperature. Thus to maintain 
the diffuse Cygnus X 2.7 GHz radio source at a distance, 

2 Qo(M) = 5.5xl0 42 M 4 s —1 , for stars born in the mass range, 
20 < M < 40 M , and Q (M) = 8.2xl0 45 M 2 s~\ 40 < M < 120 
M Q (McKee & Williams 1997). 



r, of 1.45 kpc, T e 8000 K, requires Q of 0.45 x 10 51 s _1 . 
To maintain the entire Cygnus radio emission, discrete 
sources as well as the diffuse background, requires ^10 51 
s _1 . In fact the embedded Cyg OB2 massive stars are 
potentially capable of the maintaining the entire Cygnus 
thermal radio emission. 

We propose that the ~ 6° (R ~ 75 pc) radio Cygnus 
X is the source of the N II 205 \i line emission, resolved 
by FIRAS, in the Cygnus region. Further support for 
this conjecture is provided by Lozinskaya et al. (2002). 
Approximate the 6° Cygnus X source as a sphere of ~ 3° 
radius sphere (R ^75 pc). They remark that a typical 
electron density is ^4.5 cm~ 3 found for an average line of 
sight through such a spherical region, 4R/3 ^100 pc, and 
an emission measure ~ 2000 cm~ 6 pc. Note that their 
estimated mean density is close to the mean value of 3 
cm~ 3 derived by McKee & Williams for H II envelopes. 

Thus, following the lead of Lozinskaya et al., we repre- 
sent a modeled Cygnus X source, illustrated in Figures 
4 to 6, as a uniform sphere with a 75 pc radius, whose 
center is located at 1 = 80° and b = 0° . We assume that 
the source is resolved normal to the Galactic plane. To 
approximate the effect of the 7° beam of FIRAS instru- 
ment, we weighted this source expected intensity over a 
Gaussian profile with a FWHM of 7°. To reproduce the 
observed 205 /jm intensity in Figures 4 to 6 we found 
that the required N II line flux of our modeled source is 
2.1 xl0~ 7 erg s" 1 cm~ 2 , which corresponds to a N II 205 
/im luminosity of 2.4xl0 37 erg s _1 at a distance of 1.45 
kpc. 

This powerful single source contributes 1.7x 10~ 3 of the 
Galactic 205 /im luminosity. How does this line luminos- 
ity compare to the expected Cygnus contribution to the 
total number of ionizing photons absorbed by photoion- 
ized gas in the Galaxy? Based on analyses of thermal 
radio emissions and a published study (Bennett et al. 
1994) of FIRAS N II far infrared line emissions, McKee 
& Williams estimate a total Galactic ionization rate, Qq 3 , 
of 1.9xl0 53 s _1 . However, based on their prescription 
this value should decrease by a factor, (R s /8.5) 2 ~0.8, 
due the decrease in the distance of the Sun from the 
Galactic Center of 8.5 to 7.6 kpc. However the deriva- 
tion of Q(f also relies on the value of Ljv//205/im , the total 
Galactic 205 /im luminosity. In section 5 we derived this 
luminosity from a detailed model fit to the FIRAS 205 
/im longitudinal intensities. Our 205 /Ltm luminosity es- 
timate of 1.2xl0 40 erg s~ x within the solar circle (R s 
~ 7.6 kpc) is half the published value of Bennett et al. 
(1994) of 2.4xl0 40 erg s _1 yet our smaller estimate is 
consistent with their value when ~50% uncertainties are 
included. Due to difficulties in merging disparate radio 
data sets, which affect analyses of radio observations of H 
II regions, we assign greater weight to analyses of far in- 
frared emission lines. Consequently, we estimate a total 
Galactic Q^ of 9 xlO 52 s" 1 . 

If we assume that interstellar gas in the H II shell 
around Cyg OB2 has solar composition, then its expected 
fraction of the Galactic Lyc luminosity implied by N II 
line emissions is ~ 1.78 x 1.7 x 10~ 3 = 3 x 10~ 3 . Thus 
we estimate that the value Qo required to maintain the 
Cygnus X N II line emissions to be 2.8xl0 50 s _1 . Such 
an estimated Qo is only 25% of Qo we calculated from 
the number of O stars expected to be present as well 
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as the ionizing luminosity required to maintain the total 
Cygnus X thermal emission. Consequently, we view the 
nearby Cyg OB2 as the source of N II 205 /im from the 
Cygnus region. 

7.2. Gum Nebula - Vela Region 

We interpret the bright 205 /jm FIRAS feature at 1 w 
265°, not as a spiral- arm tangent, but as a H II region 
excited primarily by a single close O star, 77-Puppis. The 
Gum Nebula - Vela Complex contains only three O stars 
compared to ^75 O stars residing in the rich Cyg OB2 
but the Complex is closer by a factor of ~ 5. Also its 
observed N II 205 fim line intensity is smaller that ^3 
than that of the Cygnus X N II source. The large, nearby 
Gum Nebula forms an Ha ring on the sky with ^18° 
radius with its center at 1 <~ 258° and b ~ -2° (Gum 
1956; Chanot & Sivan 1983). These Ha emissions most 
likely are powered by stellar UV fluxes from r\ Puppis 
(041), and 7 2 Velorum (WC8+09I binary) (e.g., Gum 
1956; Beuermann 1973; Reynolds 1976; Wallerstcin, Silk, 
& Jenkins 1980; Chanot & Sivan 1983). 

However, the origin of the Gum nebula itself is con- 
troversial. The Gum Nebula has been modeled as: 1) 
a fossil Stromgren sphere, created by the stellar progen- 
itor of the Vela supernova remnant (e.g., Brandt et al. 
1971; Alexander et al. 1971), 2) an old (>10 6 yr) super- 
nova remnant, ionized by UV fluxes from r\ Puppis and 
7 2 Velorum (Reynolds 1976), and 3) an ordinary H II re- 
gion (e.g., Gum 1956; Beuermann 1973). Subsequently, 
Wallerstein, Silk, & Jenkins (1980) investigated the prop- 
erties of interstellar gas in the Gum Nebula using both 
optical and UV line observations. They found that their 
data was generally consistent with a model of the Gum 
Nebula as an ordinary H II region, ionized by OB stars, 
rj Puppis, 7 2 Velorum, as well as by the stellar progenitor 
of the Vela supernova remnant, and the gas is stirred by 
multiple stellar winds. 

Ionized interstellar components with gas densities of 
order >1 cm~ 3 , relevant to N II 205 /im emissivities, can 
be found in a number of sites in the Gum Nebula. For 
example, Reynolds (1976) estimated that the Ha emis- 
sion shell is 15 to 30 pc thick and the electron densities 
<~ 5 - 2.5 cm -3 . Recently, Sushch, Hnatyk, and Neronov 
(2010) investigated the complex interactions of the Vela 
supernova remnant and the 5°-sized stellar wind bubble 
of 7 2 Velorum. They found that the interstellar density 
around 7 2 Velorum is ~ 12 cm~ 3 . 

The O-star catalog of Maz-Apellniz et al. (2004) lists 
the 7 2 Velorum binary at 1 = 262.8°, b = -7.7° with 
spectral types 091 and WC8 and r\ Puppis, a runaway 
O star, at 1 = 256°, b = -4.7° with spectral class 041. 
The total stellar Q is ~1.32xl0 50 ■er 1 : r\ Puppis, a very 
luminous 041 star, has an expected Qo of 1.05xl0 50 s _1 
(Vacca et al. 1996), the 091 star of 7 2 Velorum binary has 
a Qo of 2.14xl0 49 s" 1 (Vacca et al. 1996), and the WR 
star of 7 2 Velorum binary has a Qo of 6xl0 48 s _1 (Dc 
Marco et al. 2000). Using radial velocity measurements 
North et al. (2007) determined that 7 2 Velorum binary is 
close, at a distance, r, of 336±7 pc. Hipparchos parallel 
measurements provide a distance to 77 Puppis of 329^77° 
pc (van der Hucht et al. 1997). The runaway star, 77- 
Puppis contributes <~ 80% of the total stellar Q and 
these two stellar systems are close, separated only by 



~40 pc. 

Based on these parameters we assumed a single source 
of N II line emission, a uniform sphere with a 30 pc 
radius, whose center is located at 1 = 262° at a distance 
of 330 pc. To approximate the effect of the 7° beam of 
the FIRAS instrument we weighted this source intensity 
over a Gaussian profile with a full-width, half-maximum 
of 7°. We found that our model required N II line flux 
of 7.7xl0~ 8 erg s _1 cm~ 2 , which corresponds to a 205 
(Um luminosity of 10 36 erg s _1 at the assumed distance 
of 0.33 kpc. 

The Vela source contributes 10 36 /1.4x 10 40 w 7xl0~ 5 
of the total galactic 205 /im line luminosity. If we as- 
sume that interstellar gas in the vicinity of Vela has so- 
lar composition, then its expected fraction of the Galac- 
tic Lyc luminosity implied by N II line emissions is 
~ 1.78 x 7 x 10~' 5 = 1.25 x 1CT 4 . Scaling the Galactic 
Q^ by 1.25xl0 -4 , we estimate the value Qo required to 
maintain the Vela N II line emissions to be ~ 1.1 x 10 49 
s _1 . Since the Vela O-stars actually radiate ^1.3xl0 50 
s _1 , only ~10% of the stellar ionizing luminosity is re- 
quired to power the N II 205 /im emission observed by 
FIRAS. 



8. SUMMARY 

Core collapse supernovae of OB stars, formed in large 
OB associations, blow giant superbubbles, in which their 
collective shocks accelerate most of the Galactic cosmic 
rays. With the sites of cosmic ray acceleration now gen- 
erally recognized, we can now proceed to more sophis- 
ticated model tests of how the cosmic rays propagate 
within the Galaxy. In order to define the initial con- 
ditions for such models of cosmic ray propagation on a 
Galactic scale, we have constructed a three-dimensional 
spatial model of the Galactic distribution of superbub- 
bles, using the measured FIRAS Galactic longitudinal 
distribution of the 205 /im N II luminosity, emitted by H 
II envelopes surrounding hot tenuous superbubble cores. 
These H II envelopes in turn are excited by the Lyman 
continuum radiation of the massive O stars embedded in 
the superbubble cores. 

Recently Vallee (2008) critically evaluated past evi- 
dence of Galactic spiral structure and he determined up- 
to-date Galactic-plane locations of four spiral arms. We 
expanded his model by the inclusion of the fall offs of the 
spatial density of H II envelopes both parallel and nor- 
mal to the positions of spiral arms. Once suitable choices 
were made of the parameters governing these density fall 
offs we found that the observed FIRAS 205 /un N II lon- 
gitudinal distribution is reproduced well by spiral-arm 
source models once emissions from two nearby H II re- 
gions were included. One nearby source is the Cyg OB2 
association, one of the richest Galactic OB associations, 
located at a distance of 1.45 kpc, with its center at 1 = 
80.1° and b = +0.9°. The other source, at a distance 
of 0.3 kpc, is ionized by the nearest runaway O star, 
?7-Puppis and the O star binary 7 2 Velorum. Finally, 
we remark that, as illustrated by our Figures 4, 5, & 6 
the Galactic superbubble distribution, as traced by H II 
envelope emission, is dominated by the contribution of 
massive star clusters residing in the spiral arms and two 
nearby H II regions. 



10 



Higdon and Lingcnfelter 



REFERENCES 



Ackcrmann, M., Ajello, M., Allafort, A., et al. 2011, Science, 334, 
1103 

Afflcrbach, A., Churchwcll, E., & Werner, M. W. 1997, ApJ, 478, 
190 

Alexander, J. K., Brandt, J. C, Maran, S. P., et al. 1971, ApJ, 
167, 487 

Anantharamaiah, K. R. 1986, J. Astr. Ap., 7, 131 
Anderson, J. P., Habcrgham, S. M., James, P. A. et al. 2012, 

MNRAS, 424, 1372 
Azimlu, M., Marciniak, R., & Barmby, P. 2011, AJ, 142, 139 
Baade, W. 1963, in Evolution of Stars and Galaxies, (Cambridge: 

Harvard University) 
Baade, W., & Zwicky, F. 1934, Phys. Rev., 46, 76 
Barbon, R., Ciatti, F., & Rosino, L. 1979, A&A, 72, 287 
Bastian, N., Covey, K. R., & Meyer, M. R. 2010, ARA&A, 48, 339 
Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2005, ApJ, 

630, L49 

Benjamin, R. A. 2008, in Massive Star Formation: Observations 
Confronts Theory, ASP Conf. 387, eds. Beuther, H., Linz, H., 
& Henning, Th. (San Francisco: ASP) 

Bennett, C. L., Fixsen, D. J., Hinshaw, G. et al. 1994, ApJ, 434, 
587 

Beuermann, K. P. 1973, Ap&SS, 20, 27-38 
Blaauw, A. 1964, ARA&A, 2, 213 

Blaauw, A. 1991, in Physics of Star Formation and Early Stellar 

Evolution, ed. C, Lada & R. Wimmer-Schwcingruber, AIP, 257 
Bloemen, J. B., Deul, E. R., & Thaddeus, P. 1990, A&A, 233, 437 
Boulanger, F., Stark, A. A., & Combes, F. 1981, A&A, 93, LI 
Brandt, J. C, Stecher, T. P., Crawford, D. L., & Maran, S. P. 

1971, ApJ, 163, L99 
Brown, A. G. A., Hartmann, D., & Burton, W. B. 1995, A&A, 
300, 903 

Caswell, J. L.. & Haynes, R. F. 1987, A&A.171, 261 

Chanot, A., & Sivan, J. P. 1983, A&A, 121, 19 

Cohen, R. S., Dame, T. M., & Thaddeus, P. 1986, ApJS, 60, 695 

Comeron, F., & Pasquali, A. 2012, A&A, 543, A101 

Daflon, S., & Cunha, K. 2004, ApJ, 617, 1115 

Dame, T. M., Elmegreen, B. G., Cohen, R. S., et al. 1986, ApJ, 

305, 892 

de Avillez, M. A., & Breitschwerdt, D. 2005, A&A, 436, 585 
De Marco, O., et al. 2000, A&A, 358, 187 

de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J. et al. 1999, 
AJ, 117, 354 

Dessart, L., Hillier, D. J., Li, C. et si. 2012, MNRAS, 424, 2139 
Digel, S., Thaddeus, P., & Bally J. 1990, ApJ, 357, L29 
Dobbs, C. L., Bonncll, I. A., & Pringle, J. E. 2006, MNRAS, 
371,1663 

Downes, D., Wilson, T. L., Bicging, J., et al. 1980, A&AS, 40, 379 
Drilling, J. S. & Landolt, A. U. 2000, p. 389 in Allen's 

Astrophysical Quantities, ed. A. N. Cox (4th ed; NY: Springer) 
Drimmel, R. & Spergel, D. N. 2001, ApJ, 556, 181 
Drew, J. E., Dreimel, R., Irwin, M. J. et al. 2008, MNRAS, 386, 

1761 

Eferemov, Y. N. 2011, Astronomy Rcoprts, 55, 108 
Englmaier, P., & Gerhard, O. 1999, MNRAS, 304, 512 
Filippenko, A. V. 1997, ARA&A, 35, 309 

Fixsen, D. J., Bennett, C. L., & Mather, J. C. 1999, ApJ, 526, 207 
Forbes, D. 1983, in Kinematics, Dynamics and Structure of the 

Milky Way, ed. W. H. L. Shutcr (Ricdcl: Dordrecht, NL), 

217-22 

Gaskell, C. M., Cappellaro, E., Dinerstein, H. L. ct al. 1986, ApJ, 

306, L77 

Georgelin, Y. M., & Georgelin, Y. P. 1976, A&A, 49, 57 
Gies, D. R. 1987, ApJS, 64, 545 

Ginzburg, V. L., & Syrovatskii, S. I. 1964, The Origin of Cosmic 

Rays (New York: Macmillan) 
Georgy, C, Meynet, G., Walder, R. et al. 2009, A&A, 502, 611 
Grabelsky, D. A., Cohen. R. S., Bronfman, L., ct al. 1988, ApJ, 

331, 181 

Gum, C. S. 1956, Observatory, 17, 150 
Hanson, M. M. 2003, ApJ, 597, 957 
Hayakawa, S. 1956, Prog. Theor. Phys., 15, 111 
Hayakawa, S., Matsumoto, T., Murakami, H., et al. 1981, A&A, 
100, 116 



Heiles, C, Haffner, L. M., & Reynolds, R. J. 1999 in ASP Conf. 
Ser. 168 New Perspectives on the Interstellar Medium, ed. A. 
R. Taylor, T. L. Landcckcr, & G. Jonas (San Francisco, CA: 
ASP), 211 

Heiles, C. 2001 in ASP Conf. Ser. 231 Tetons 4: Galactic 

Structure, Stars, & Interstellar Medium, ed. C. E. Woodward, 
M. D. Biscay, & J. M. Shull (San Francisco, CA: ASP), 294 
Heyer, M. H., Carpenter, J. M., & Snell, R. L. 2001, ApJ, 551, 852 
Higdon, J. C, & Lingenfelter, R. E. 2003, ApJ, 590, 822 
Higdon, J. C, & Lingenfelter, R. E. 2005, ApJ, 628, 738 
Higdon, J. C, Lingenfelter, R. E., & Ramaty, R. 1998, ApJ, 509, 
L33 

Hou, L. G., Han, J. L., & Shi, W. B. 2009, A&A, 499, 473 

Huchtmeier, W. K., & Wendker, H. J. 1977, A&A, 58, 197 

Humphreys, R. M. 1978, ApJS, 38, 309 

Ikhsanov, R. N. 1961, SvA, 4, 923 

James, P. A., & Anderson, J. P. 2006, A&A, 453, 57 

Kapteyn, J. C. 1914, ApJ, 40, 43 

Kennicutt, R. C, & Hodge, P. W. 1976, ApJ, 207, 36 

Kiminki, D. C, Kobulnicky, A., Kincmuchi, K. ct al. 2007, AJ, 

664, 1102 
Knodlseder, J. 2000, A&A, 360, 539 

Korpi, M. J., Brandenburg, A., Shukurov, A., ct al. 1999, A&A, 
350, 230 

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 
Lee, J. H., Hwang, N., & Lee, M. G. 2011, ApJ, 735, 75 
Lcitherer, C, Schaerer, D., Goldader, J. D. et a.l, 1999, ApJS, 
123, 3 

Levine, E.S., Blitz, L., & Heiles, C. 2006, Science, 312, 1773 
Lozinskaya, T. A., Pravdikova, V. V, & Finoguenov, A. V. 2002, 
AstL, 28, 223 

Mai'iz-Appellaniz, J., Walborn, N. R., Galuc, H. A., ct al. 2004, 
ApJS, 151, 103 

Martin, D. C, Fanson, J., Schiminovich, D., et al. 2005, ApJ, 
619, LI 

Mason, B. D., Gies, D. R., Hartkopf, W. I. 1998, AJ, 115, 821 
Massey, P. 2003, ARA&A, 41, 15 

Masscy, P., & Meyer, M. R. 2001, Encyclopedia of Astronomy 

and Astrophysics, (NY: Nature Publishing Group) 
Massey, P., & Thompson, A. B. 1991, AJ, 101, 1408 
May, J., Alvarez, H., & Brofman, L. 1997, A&A, 327, 325 
McCray, R. & Snow, T. P. 1979, ARA&A, 17, 213 
McCray, R. & Kafatos, M. 1987, ApJ, 317, 190 
McKec, C. F. & Willams, J. P. 1997, ApJ, 476, 144 
Minkowski, R. 1941, PASP, 53, 224 

Moffat, A. F. J., Fitzgerald, M. P., & Jackson, P. D. 1979, ApJS, 
38, 197 

Morgan, W. W., Sharpless, S., & Osterbrock, D. 1952, AJ, 57, 3 
Nakagawa, T., Yui, Y. Y., Doi, Y. ct al. 1998, ApJS, 115, 259 
Nakanishi, H., & Sofue, Y. 2003, PASJ, 55, 191 
Nakanishi, H., & Sofue, Y. 2006, PASJ, 58, 847 
Ncugent, K. F., & Massey, P. 2011, ApJ, 733, 123 
North, J. R., Tuthill. P. G., Tango, W. J., et al. 2007, MNRAS, 
377, 415 

Ogliore, R. C, Stone, E. C, Leske, R. A., et al. 2009, ApJ, 695, 
666 

Paladini, R., Davies, R. D., & DeZotti, G. 2004, MNRAS, 347, 
237 

Paczyhki, B. 1967, Acta Astronomica, 17, 355 
Podsiadlowski, Ph., Joss, P. C, & Hsu, J. J. L. 1992, ApJ, 391, 
246 

Rauch, B. F., Link, J. T., Lodders, K., ct al. 2009, ApJ, 697, 2083 

Reed, B. C. 2000, AJ, 120, 314 

Reynolds, R. J. 1976, ApJ, 203, 151 

Reynolds, R. J., & Ogden, P. M. 1979, ApJ, 229, 942 

Roshi, D. A., & Anantharamaiah, K. R. 2000, ApJ, 557, 226 

Rubin, R. H. 1985, ApJS, 57, 349 

Rumstay, K. S., & Kaufman, M. 1983, ApJ, 274, 611 

Russeil, D. 2003, A&A, 397, 133 

Russeil, D. Adami, C, Georgelin, Y. M. 2007, A&A, 470, 161 
Salpeter, E. E. 1955, ApJ, 121, 161 

Shaver, P. A., McGec, R. X., Newton, L. M., Danks, ct al. 1983, 

MNRAS, 204, 53 
Smartt, S. J. 2009, ARA&A, 47, 63 

Smartt, S. J., Eldridge, J. J., Crockett, R. M., et al. 2009, 
MNRAS, 395, 1409 



Galactic Distribution Supcrbubbles 



11 



Smith, L. F, Biermann, P., & Mezger, P. G. 1978, A&A, 66, 65 
Smith, N., Li, W., Filippenko, A. V., et al. 2011, MNRAS, 412, 
1522 

Solomon, P. M., Sanders, D. B., & Rivolo, A. R. 1985, ApJ, 292, 
L19 

Stark, A. A., & Lee, Y. 2005, ApJ, 619, L159 
Stark, A. A., & Lee, Y. 2006, ApJ, 641, L113 
Steiman-Cameron, T. Y., Wolfire, M., & Hollcnbach, D. J. 2010, 

ApJ, 722, 1460 
Stil, J., Wityk, N., Ouyed, R., ct al. 2009, ApJ, 701, 330 
Sushch, I., Hnatyk, B., & Neronov, A. 2010, A&A, 525, A154 
Taylor, J. H., & Cordes, J. M. 1993, ApJ, 411, 674 
Tenorio-Tagle, G. & Bodenheimer, P. 1988, ARA&A, 26, 145 
Tomisaka, K., & Ikeuchi, S. 1986, PAS J, 38, 687 
Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914 
Vallee, J. P. 2008, AJ, 135, 1301 

Van der Hucht, K. A. et al. 1997, New Astronomy, 2, 245 



Vazquez, R. A., May, J., Bronfman, L., ct al. 2008, ApJ, 672, 930 
Veron, P. 1965, Annalcs d.Astrophysiquc, 28, 391 
Wallerstein, G., Silk, J., & Jenkins, E. B. 1980, ApJ, 240, 834 
Watson, C, Araya, E., Sewilo, M. et al. 2003, ApJ, 587, 714 
Weaver, H. 1970 in IAU 38, Spriral Structure of our Galaxy, ed. 

W. Becker & G. I. Kontopoulos (Dordrecht: Reidel), 126 
Wendker, H. J. 1970, A&A, 4, 378 
Williams, J. P., & McKee, C. F. 1997, ApJ, 476, 166 
Wright, E. L., Mather, J. C, Bennett, C. L, et al. 1991, ApJ, 381, 

200 

Wright, N. J., Drake, J. J., Drew, J. E., & Vink, J. S. 2010, ApJ, 
713, 871 

Zatsepin, V. I., & Sokolskaya, N. V. 2007, AstL, 33, 25 
Zwicky, F., 1939, Phys. Rev., 55, 986 



12 Higdon and Lingenfelter 



i.5xicr 4 



1.0x10" 



2 5.0xl0" : 



0.0 



ill 



1 



ri, 



U p = 2.5 kpc & 3 < p < 10 kpc 
Axisymmetric Disk 




A 




180 150 120 90 60 30 330 300 270 240 210 180 
Galactic Longitude 



Figure 1. Axisymmetric Disk Source Model. A modeled N II 205 fim longitudinal intensity distribution based on a radially truncated 
axisymmetric disk population with H p = 2.5 kpc and R s = 7.6 kpc. A Gaussian distribution transverse the Galactic plane with dispersion, 
a z , of 0.15 kpc is assumed. The modeled intensity distribution is normalized to the observed FIRAS intensity at 1 = 30° requiring a total 
Galactic N II 205 (im luminosity of 1.85xl0 40 erg s — 1 . Shown as a histogram are the mean longitudinal FIRAS 205 (im intensities, Fig. 
5e (Fixsen, Bennett, & Mather 1999) averaged over Galactic longitudes of Al = 5° and over Galactic latitudes (centered on the Galactic 
plane) of Ab = 1°. 
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Figure 2. Our model of the Galactic distribution of the median of spiral arms based on the prescriptions of Valllee (2008). The thick 
solid curves represent equation (1), for Perseus, Norma-Cygnus, Crux-Scutum, and Sagittarius-Carina arms. Mean longitudes as well as 
root-mean-square uncertainties of measured spiral-arm tangents from Table 2 of Vallee (2008) are represented by thin solid lines. The 
curves to either side of the arms represent the one cr Gaussian 0.5 kpc half-arm widths. 
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Figure 3. Representative 205 fim N II intensities calculated from Equation (10) for a total Galactic spiral-arm luminosity of Ljv II 205 um 
of 1.4xl0 40 erg s — 1 and nominal arm parameters of H p = 2.4 kpc and a a = 0.5 kpc are displayed. Four different sets of relative arm 
luminosities, fg^, fgp, fjvCi an d fp, have been considered and we find that the second set of relative luminosities appears to give the best 
fit of the four, so as we discuss in the text, we use them hereafter. As in Fig. 1 the mean longitudinal FIRAS 205 fim intensities, Fig. 5c 
(Fixscn, Bennett, & Mather 1999) averaged over Galactic longitudes of Al = 5° are shown as a histogram. 
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Figure 4. Longitudinal Spiral-Arm Distribution as Function of Variations in H p . Modeled representative 205 fira N II intensities calculated 
for a total spiral-arm luminosity of Ljy n 205 pm °f 1.4xl0 40 erg s — 1 and a fixed arm-width parameter, a a of 0.5 kpc are plotted for 
a range of radial parameters, H p , between 1.8 to 3.0 kpc, showing the value of about 2.4 appears to give the best fit and we use that 
hereafter. As in Fig. 1 the mean longitudinal FIRAS 205 fim intensities, Fig. 5e (Fixsen, Bennett, &; Mather 1999) averaged over Galactic 
longitudes of Al = 5° are shown as a histogram. Also included, in addition to the modeled spiral arms, are the contributions of nearby 
local H II regions, Cygnus X, at 1 ?s 80° and Vela Complex, at 1 262°. 
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Figure 5. Longitudinal Spiral-Arm Distribution as Function of Variations in a a- Modeled representative 205 £tm N II intensities calculated 
for a Ljv n 205 um °f 1.4xl0 40 erg s _1 and a fixed radial parameter, H p of 2.4 kpc for a range of arm-width parameters, a a, between 0.25 
to 0.75 kpc. We find that a Gaussian half-width of 0.5 ±0.1 kpc appears to give the best fit and we use that hereafter. As in Fig. 1 the 
mean longitudinal FIRAS 205 ^m intensities, Fig. 5e (Fixsen, Bennett, & Mather 1999) averaged over Galactic longitudes of Al = 5° are 
shown as a histogram. Also included, in addition to the modeled spiral arms, are the contributions of nearby local H II regions, Cygnus X, 
at 1 m 80° and Vela Complex, at 1 262° . 
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Figure 6. Individual Spiral-Arm Contributions to 205 fim N II Intensities. The solid curve shows the modeled 205 fim N II intensity 
distribution for a total spiral-arm luminosity, Ljv II 205 umi °f 1.4xl0 40 erg s — 1 and a fixed arm-width parameter, a a, of 0.5 kpc for a 
radial parameter, H p , of 2.4 kpc. As in Fig. 1 the mean longitudinal FIRAS 205 fim intensities, Fig. 5e (Fixsen, Bennett, & Mather 1999) 
averaged over Galactic longitudes of Al = 5° are shown as a histogram. Also included, in addition to the modeled spiral arms, are the 
contributions of nearby local OBA, Cygnus X, at 1 ss 80° and Vela Complex, at 1 Ri 262° . 
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Figure 7. A representative planar spiral-arm Lyc source distribution which we also expect to reflect the Galactic spatial distribution of 
OB associations and their surrounding superbubbles. A Galactic-centered Cartesian grid with spacings dx = 0.2 kpc and dy = 0.2 kpc was 
used to create the illustrated arm distribution, corresponding to a nominal H p = 2.4 kpc and a a = 0.5 kpc. It was calculated employing 
the Sx (Pi 6) interpolants, discussed in the text, and the simple transformation x = pcos6> and y =psin0. Note that fS(p, 9)pdpd9 = 1. For 
purposes of illustration, the solar position is identified by a white square. 



